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SUMMARY 


Shockwave  and  turbulent  boundary  layer  interactions  produce  intense  localized  pressure  loads 
and  heating  rates  that  can  have  a  dramatic  influence  on  the  drag  and  heating  experienced  by  a 
high-speed  vehicle,  and  can  significantly  impact  fuel  mixing  and  combustion  in  propulsion 
systems.  The  lack  of  standardized  and  traceable  databases  prevents  the  calibration  of 
computational  fluid  dynamic  models  to  accurately  represent  these  critical  flow  phenomena. 

The  effort  described  here  is  the  result  of  a  one-year  research  program  FA9550-09- 1-0464  and 
provides  accurate  statistical  evidence  of  the  upstream  and  downstream  influence  in  shock 
unsteadiness  in  a  fully  separated  configuration  of  boundary  layer  flow  over  a  compression  comer 
at  Mach  2.9  and  Re0  =2900.  Direct  numerical  simulation  data  was  gathered  over  1000  S/Ux 
(350  Lsep  /[/„),  making  the  analysis  of  the  low-frequency  unsteadiness  possible.  The  dominant 

scales  in  the  flow  were  investigated  by  spectral  analysis,  and  the  statistical  link  between  the  low- 
frequency  shock  motion  and  the  upstream  and  downstream  flow  were  investigated.  Strong 
coherence,  at  the  low  frequencies  of  the  shock  motion,  is  observed  with  the  downstream, 
separated  flow.  The  coherence  with  the  upstream,  undisturbed  boundary  layer  is  shown  to  be 
statistically  significant,  but  weak.  The  findings  were  reported  in 

■  Priebe,  S.,  and  Martin,  M.P.,  “Direct  Numerical  Simulation  of  Shockwave  and  Turbulent  Boundary  Layer 
Interactions”,  AIAA  Paper  2009-0589 ,  47th  AIAA  Aerospace  Science  Meeting  and  Exhibit ,  Orlando,  FL, 
January  2009. 

A  direct  numerical  simulation  of  a  reflected  shock  and  turbulent  boundary  layer  interaction  was 
also  characterized  using  DNS  during  this  period  effort  and  appears  in: 

■  Priebe,  S.,  Wu,  M.,  and  Martin,  M.P.,  “Direct  Numerical  Simulation  of  a  Reflected  Shockwave  and 
Turbulent  Boundary  Layer  Interactions”,  AIAA  Journal,  1,  5,  2009. 

The  above  work  builds  upon  AFOSR  previously  sponsored  research  FA9550-090- 1-0464,  where 
we  accomplished  the  development  and  validation  against  experiments  at  the  same  flow  and 
boundary  conditions  of  direct  numerical  simulations  of  shock  and  turbulent  boundary  layer 
interactions.  We  pioneered  the  development  of  a  unique  numerical  capability  that  allows  the 
accurate  and  detailed  three-dimensional  turbulence  data  at  a  reasonable  tum-around  time.  In 
turn,  parametric  studies  of  fundamental  flow  physics  are  feasible,  for  the  first  time.  By  accurate, 
it  is  meant  that  the  numerical  uncertainty  is  within  the  experimental  error.  By  reasonable  tum- 
around  time,  it  is  meant  that  the  computational  time  is  comparable  to  the  experimental  turn¬ 
around-time.  The  numerical  methods,  the  simulations  and  their  validation  against  experimental 
data  have  been  published  in  the  following  journal  papers: 

•  Taylor,  E.M.,  Wu,  M.,  and  Martin,  M.P.,  “Optimization  of  Nonlinear  Error  Sources  for  Weighted  Non- 
Oscillatory  Methods  in  Direct  Numerical  Simulations  of  Compressible  Turbulence,”  Journal  of 
Computational  Physics ,  223,  384-397,  2007. 

•  Wu,  M.,  and  Martin,  M.P.,  “Direct  Numerical  Simulation  of  Shockwave  and  Turbulent  Boundary  Layer 
Interaction  induced  by  a  Compression  Ramp,”  AIAA  Journal ,  45,  4,  879-889,  2007. 

•  Ringuette,  M.,  Wu,  M.,  and  Martin,  M.P.,  “Coherent  Structures  in  DNS  of  Turbulent  Boundary  Layers  at 
Mach  3  f  Journal  of  Fluid  Mechanics,  594,  59-69,  2008 
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The  unsteady  motion  of  STBLI  has  been  analyzed  using  the  DNS  data  and  this  work  has  been 
published  in: 

•  Wu,  M.  and  Martin,  M.P.,  “Analysis  of  Shock  Motion  in  STBLI  using  Direct  Numerical  Simulation  Data,” 
Journal  of  Fluid  Mechanics,  594,  71-83,  2008. 

In  addition,  the  data  analysis  using  the  DNS  of  Wu  &  Martin  suggest  that  low-Reynolds  number 
shock-wave  turbulent  boundary  layer  interactions  exhibit  differences  with  previous 
measurements  at  high  Reynolds  number.  The  low  Reynolds  number  effects  are  due  to  the  greater 
influence  of  viscosity,  and  result  in  a  smaller  peak  in  the  RMS  of  the  wall  pressure  fluctuations, 
an  enriched  intermittency  of  the  wall-pressure  signal,  and  a  substantially  larger  separation  zone. 
Unlike  previous  studies  at  high  Reynolds  number,  the  richer  wall-pressure  signal  of  the  low- 
Reynolds  number  data  cannot  be  used  to  determine  the  location  of  the  shock  wave.  The  primary 
shock  wave  does  not  penetrate  as  deeply  into  the  boundary  layer  as  for  the  high  Reynolds 
number  flows,  so  it  is  more  accurate  to  determine  the  low-Reynolds  number  shock  location  in  the 
outer  region  of  the  boundary  layer.  Despite  the  difference  between  the  low  and  high  Reynolds 
number  data,  the  low-frequency  shock  motion  (relative  to  the  high-frequency  that  characterizes 
the  undisturbed  boundary  layer)  reported  for  high  Reynolds  number  flows,  and  the  turbulence 
amplification  across  the  interaction  region,  are  not  affected  by  the  low  Reynolds  number 
condition.  These  findings  have  been  published  in  the  following  journal  articles: 

•  Wu,  M.,  and  Martin,  M.P.,  “Direct  Numerical  Simulation  of  Shockwave  and  Turbulent  Boundary  Layer 
Interaction  induced  by  a  Compression  Ramp,”  AIAA  Journal ,  45,  4,  879-889,  2007. 

•  Ringuette,  M.J.,  Wu,  M.,  and  Martin,  M.P.,  “Low  Reynolds  Number  Effects  in  a  Mach  3  Shock  Turbulent 
Boundary  Layer  interaction,”  AIAA  Journal ,  46,  7,  2008. 

Also,  as  part  of  this  effort,  accurate  wall-pressure  data  for  Mach  3  interactions  at  low  Reynolds 
numbers,  accessible  to  DNS  and  LES,  have  been  gathered  and  published  in  the  following  journal 
article: 

•  Ringuette,  M.J.,  Bookey,  P.,  Wychham,  C.,  Smits,  A.J.,  “Experimental  Study  of  a  Mach  3  Compression 
Ramp  Interaction  at  Re0  =  2400”,  AIAA  Journal ,  47,  2  2009. 
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1.  Background 


1.1.  Canonical  configurations 

The  simplest  canonical  shock  and  turbulent  boundary  layer  configurations  are  shown  in  Fig.  6, 
namely  a  compression  comer,  a  reflected  shock  interaction,  and  a  sharp  fin  interaction.  The  flow 
features  for  these  configurations  and  our  learning  from  studying  these  interactions  are  briefly 
summarized  below. 


Compression  corner  interaction 


The  compression  comer  interaction  is  one  of  the  simplest  cases  of  STBLI  that  occur  in  internal 
and  external  vehicle  flows.  This  configuration  has  been  extensively  studied  experimentally  by, 
for  example,  Settles  et  al  (1979),  Kuntz  et  al.  (1987),  Smits  &  Muck  (1987),  Dolling  and  Murphy 
(1983),  Ardonceau  (1984),  and  Selig  et  al.  (1989).  The  early  research  covers  a  wide  range  of 
turning  angles  and  Reynolds  numbers,  where  the  lowest  Reynolds  number  reported  is  Re@  = 
23,000  (Settles  et  al,  1978),  and  the  upper  Mach  number  is  limited  to  about  5,  corresponding  to 
the  experiments  by  Erengil  &  Dolling  (1991). 

We  have  learned  much  from  the  high-Reynolds  number  experiments.  The  pressure  gradient 
imposed  by  the  shock  can  cause  the  flow  to  separate  in  the  vicinity  of  the  comer  location,  and  at 
Mach  2.9  the  flow  is  on  the  verge  of  separation  with  a  comer  angle  of  16°  (called  incipient 
separation).  At  24°,  the  time-averaged  region  of  separation  spans  about  25,  starting 
approximately  1.2  5  ahead  of  the  comer  and  reattaching  at  about  0.8  5  downstream  of  the  comer. 
Near  the  line  of  separation,  compression  waves  merge  into  a  well-defined  separation  shock,  and 
a  second  shock  forms  near  the  line  of  attachment  (Settles,  1976).  Figure  7  illustrates  the  shock 
system  in  a  compression  comer  configuration  with  increasing  compression  angle.  The 
corresponding  wall  pressure  distribution  shows  an  inflection  point  or  “plateau”  in  the  region  of 
separation,  as  shown  in  Figure  8.  Further  downstream,  the  wall-pressure  eventually  recovers  to 
the  inviscid  oblique-shock  value,  but  the  point  where  this  occurs  is  located  farther  downstream 
with  increasing  compression  angle  (Settles  et  al.  1978).  For  the  24°  case,  the  inviscid  value  is 
not  recovered  before  the  end  of  the  experimental  model  is  reached,  nearly  eight  boundary  layer 
thicknesses  downstream  of  the  comer. 

A  measure  of  the  upstream  influence  is  the  distance  from  the  comer  at  which  the  shock  presence 
is  first  felt.  A  measure  of  the  streamwise  interaction  extent  is  the  separation  length,  being  the 
distance  between  the  separation  and  reattachment  points.  These  characteristic  lengths  are 
determined  from  time-averaged  measurements,  and  they  vary  with  time  due  to  the  highly 
unsteady  motion  of  the  separation  shock.  The  distance  over  which  the  separation  shock  moves 
increases  with  turning  angle,  and  at  24°  it  moves  about  0.55  (Selig  et  al.,  1989).  The  frequency  is 
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typically  an  order  of  magnitude  lower  than  any  characteristic  turbulence  frequencies.  Thus,  the 
frequency  and  scale  of  the  shock  motion  are  needed  to  fully  characterize  the  interaction. 

The  influence  of  the  compression  on  the  turbulence  is  an  enhanced  mixing  due  to  the  formation 
of  large-scale  eddies  (Kuntz  et  al,  1987)  as  the  incoming  boundary  layer  is  driven  out  of 
equilibrium.  The  boundary  layer  mean  flow  recovery  distance  increases  with  increasing 
interaction  strength  (Smits  and  Muck,  1987;  Selig  et  al.,  1989;  Ardonceau,  1984;  Kuntz  et  al, 
1987).  The  turbulence  levels  are  strongly  amplified  across  the  shock  system,  and  Selig  et  al., 
(1989)  found  that  at  Mach  2.9  the  mass-flux  fluctuations  increased  by  more  than  a  factor  of  four 
with  a  24°  turning  angle.  The  flow  distortion  is  also  seen  in  the  heat  transfer:  Evans  &  Smits 
(1996)  found  that  the  Reynolds  analogy  factor  increased  by  a  factor  of  three  through  a  16° 
interaction,  and  showed  little  sign  of  relaxation  downstream  of  the  comer. 

In  contrast  with  numerous  experimental  data,  there  are  few  detailed  simulations  such  as  DNS 
(Adams,  2000;  Martin  &  Wu,  2007)  and  LES  and  hybrid  LES/RANS  (Rizzetta  &  Visbal,  2001; 
Loginov  et  al.,  2006;  Edwards,  Choi  &  Boles,  2008).  Recently,  a  number  of  experiments 
(Bookey  et  al.,  2005;  Ringuette  &  Smits,  2007)  have  been  performed  at  lower  DNS-  and  LES- 
accessible  Reynolds  numbers,  so  that  the  validation  of  DNS  and  LES  is  possible.  From  the 
validated  DNS  data  at  low- Reynolds  numbers  we  have  learned  that  greater  viscous  effects 
diffuse  the  shock  near  the  wall,  which  results  in  a  reduced  magnitude  of  the  peak  in  the  RMS  of 
the  wall-pressure  fluctuations,  a  richer  intermittency  of  the  wall-pressure  signal,  and  a  greater 
spread  of  the  separation  zone  (Wu  &  Martin,  2007;  Ringuette,  Wu  &  Martin,  2008).  While  in 
high  Reynolds  number  experiments  the  shock  motion  is  inferred  from  measurements  of  the  wall- 
pressure,  for  low  Reynolds  number  flows  the  single  shock  does  not  penetrate  as  deeply  into  the 
boundary  layer  and  the  shock  location  is  not  well  defined  in  the  lower  half  of  the  boundary  layer. 
Thus,  at  low  Reynolds  numbers,  the  shock  motion  is  more  accurately  studied  in  the  outer  part  of 
the  boundary  layer  and  in  the  ffeestream  (Wu  &  Martin,  2007).  Low-Reynolds  number  effects, 
however,  do  not  alter  the  characteristic  low-frequency  unsteadiness  of  the  shock  wave  and  the 
separation  bubble  or  the  magnitude  of  the  turbulence  amplification  across  the  interaction  (Wu  & 
Martin,  2008;  Ringuette,  Wu  &  Martin,  2008). 

Experimental  data  at  low  Reynolds  number  and  Mach  8  for  an  8°  compression  comer  flow  have 
been  recently  reported  (Bookey  et  al.,  2005).  There  is  much  to  be  learned  from  this  flow,  even 
regarding  the  validity  of  the  strong  Reynolds  Analogies  in  the  Mach  8  turbulent  boundary  layer 
incoming  for  this  interaction  configuration. 


Re/7ecteJ  shock  interaction 


This  type  of  interaction  has  not  been  studied  as  extensively  as  the  compression  comer  case.  A 
review  of  the  reflected  shock  interaction  is  given  in  Delery  and  Marvin  (1986).  The  strength  of 
the  incident  shock  determines  the  nature  of  the  interaction.  For  weak  incident  shocks,  as 
illustrated  in  Figure  9,  the  result  is  close  to  the  inviscid  interaction.  The  incident  shock  (Cl) 
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curves  progressively  as  it  penetrates  the  boundary  layer  due  to  the  decrease  in  Mach  number 
within  the  layer.  The  effect  of  the  incident  shock  is  felt  upstream  by  pressure  propagation 
through  the  subsonic  part  of  the  boundary  layer  very  near  the  wall.  The  thickened  subsonic 
region  generates  outgoing  compression  waves  that  coalesce  into  the  reflected  shock  (C2). 

When  the  incident  shock  is  strong  enough,  the  boundary  layer  separates  and  the  flow  no  longer 
resembles  the  inviscid  case  (see  Figure  10).  Here,  the  boundary  layer  separates  at  point  S. 
Compression  waves  emanate  from  the  separation  region  and  coalesce  into  the  separation  shock 
(C2).  The  separation  shock  intersects  the  incident  shock  at  H  and  generates  the  refracted  shocks 
C3  and  C4.  Shock  C3  enters  the  boundary  layer  and  reflects  off  the  separated  region  into  an 
expansion  fan.  The  fan  turns  the  flow  towards  the  wall,  which  decreases  the  height  of  the 
separation  bubble  until  the  flow  reattaches  at  point  R.  The  resulting  reattachment  compression 
waves  compress  the  flow  gradually. 

Green  (1970)  observed  the  similarities  of  the  compression  comer  and  reflected  wave 
configurations  and  pointed  out  that  a  compression  comer  of  angle  20  will  produce  the  same 
series  of  compression  interactions  at  separation  and  reattachment  as  an  incident  shock 
configuration  with  initial  deflection  angle  0.  In  turn,  the  overall  pressure  change  is  the  same  and 
the  surface  pressure  distributions  are  nearly  identical  (Shang  et  al,  1976).  The  scaling  of  the 
upstream  influence  length  and  separation  length  is  expected  to  behave  similarly  to  the 
compression  comer  interaction.  One  major  difference  between  the  compression  comer  and  the 
incident  shock  cases  is  a  larger  separation  bubble  for  the  incident  shock  configuration,  and,  in 
contrast  to  a  compression  comer  flow  where  the  separation  bubble  height  is  only  a  small  fraction 
of  the  incoming  boundary  layer  thickness,  the  height  of  the  separation  region  in  a  reflected  wave 
flow  is  typically  comparable  to  the  incoming  boundary  layer  thickness.  Thus,  the  velocity 
profiles  show  significant  regions  of  flow  reversal,  which  result  in  a  significantly  increased 
boundary  layer  thickness  over  the  separated  region  (Delery  and  Marvin,  1986). 

More  recently,  Dupont,  Haddad  and  Debieve  (2006)  and  Dussauge,  Dupont  and  Debieve  (2006) 
studied  a  reflected  shock  interaction  at  M  =  2.3  and  Ree  =6900  generated  by  an  oblique  shock 
wave  with  deflection  angles  varying  from  7°  to  9.5°.  Detailed  Mach  number,  velocity,  and 
fluctuation  profiles  were  obtained.  The  principal  observations  agreed  with  earlier  studies  on 
separation  length  and  recovery  distance,  and  significant  unsteadiness  was  observed.  The  three- 
dimensional  nature  of  the  interaction  was  also  noted,  with  strong  swirling  motions  appearing  to 
terminate  the  separation  zone  in  the  spanwise  direction.  Bookey  et  al.  (2005)  report 
experimental  data  for  a  12°  incident  shock  on  a  boundary  layer  at  Mach  2.9  and  Ree=2400, 
providing  surface  pressure  distributions  and  Pitot  surveys,  and  flow  visualization  using  surface 
oil  and  CO2  enhanced  filtered  Rayleigh  scattering  (FRS). 

Again,  there  are  few  detailed  simulations  such  as  DNS  (Pirozzoli  &  Grasso,  2006;  Priebe,  Wu,  & 
Martin  2008)  and  LES  (Gamier,  Sagaut,  and  Deville,  2002).  Direct  comparison  against 
experiments  for  this  configuration  is  not  possible,  where  significant  three-dimensional  effects 
due  to  wind  tunnel  walls  are  observed  experimentally  (Bookey  et  al.,  2005;  Dussauge,  Dupont,  & 
Debieve,  2006).  Figure  1 1  shows  surface  oil  visualizations  (Bookey  et  al.  2005),  illustrating  the 
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significant  spanwise  variation  of  the  flow.  DNS  including  the  experimental  spanwise  length, 
wind  tunnel  side  wall  boundary  conditions  are  not  feasible  and  LES  have  not  yet  been  validated 
for  these  type  of  flows.  Robust  DNS  methodologies  for  STBLI,  however,  have  been  validated 
for  the  compression  comer  interaction  against  experiments.  Using  such  methods,  Priebe,  Wu  & 
Martin  (2009)  report  the  detailed  DNS  data  on  an  incident  shock  configuration  that  is  generated 
by  a  12°  wedge  in  the  free  stream  of  a  Mach  2.9,  Ree=2390  turbulent  boundary  layer.  They 
report  on  the  evolution  of  the  mean  and  fluctuating  flow  quantities,  the  validity  of  the  Strong 
Reynolds  analogies  and  the  characteristic  low  frequency  of  shock  motion. 


Sharp  fin  irtduced_  interaction 

A  swept-shock  interaction  is  generated  when  a  sharp  fin  placed  at  an  angle  of  attack  to  the 
incoming  flow.  Here  the  oblique  shock  sweeps  across  the  incoming  boundary  layer,  and  strong 
secondary  flows  can  be  produced  by  the  spanwise  pressure  gradients.  As  the  boundary  layer 
enters  the  rising  pressure,  the  gradient  in  Mach  number  inside  the  boundary  layer  will  cause  the 
flow  near  the  wall  to  turn  through  a  greater  angle  than  the  flow  away  from  the  wall  (as  long  as 
the  pressure  gradient  dominates).  The  differential  turning  leads  to  a  helical  secondary  flow. 
Typically,  one  or  more  large-scale  vortical  motions  are  induced  which  sweep  the  low-momentum 
fluid  from  the  near-wall  region  of  the  incoming  boundary  layer  in  the  direction  along  the  shock 
(see  Figure  12).  The  high  momentum  fluid  in  the  outer  part  of  the  boundary  layer  passes  over 
the  vortex  with  a  turning  angle  more  typical  of  the  inviscid  deflection  associated  with  the  shock, 
and  it  is  then  swept  close  to  the  wall.  The  skin  friction  and  heat  transfer  levels  seem  largely 
unaffected  by  the  strong  secondary  motions,  but  the  values  rise  sharply  in  the  region  closer  to  the 
fin  where  the  high  momentum  fluid  “scours”  the  wall.  The  turbulence  response  is  not  well 
understood.  Very  few  experimental  results  are  available,  but  measurements  by  Tan  (see  Smits  & 
Dussauge,  2006)  suggest  that  the  turbulence  levels  are  strongly  amplified,  and  Tran  et  al.  (1985) 
found  that  the  shock  is  unsteady,  leading  to  strong  wall-pressure  fluctuations.  In  these  respects, 
three-dimensional  interactions  appear  to  be  similar  to  their  two-dimensional  counterparts,  but  the 
detailed  response  of  the  turbulence  is  quite  different.  In  particular,  the  turbulence  amplification 
and  the  unsteady  pressure  loading  are  weaker. 

The  shock  bifurcates  in  response  to  the  formation  of  the  helical  secondary  flow,  in  a  manner  very 
similar  to  that  seen  in  two-dimensional  separated  compression-comer  flows.  There  is  an  initial 
turning  and  compression  by  a  well-defined  shock,  which  is  slanted  forward  (the  “separation” 
shock),  and  a  stronger  trailing  shock,  where  the  two  shock  structures  encompass  the  large-scale 
vortical  flow.  When  the  flow  is  viewed  along  the  axis  of  the  helix  it  appears  similar  to  the  cross- 
section  of  a  two-dimensional  separated  flow.  In  that  view,  a  bubble-type  separated  flow  is 
observed,  and  the  flow  characteristics  typically  scale  in  conical  coordinates.  The  experiments 
show  that  the  wall  pressure  distribution  and  the  total  pressure  distribution  can  be  collapsed  in 
conical  coordinates.  One  feature  that  deserves  particular  attention  is  the  “impinging  jet”,  found 
in  close  proximity  to  the  fin  itself.  As  the  model  by  Garg  and  Settles  (1993)  makes  clear,  the  jet 
is  formed  by  high-momentum  fluid  from  the  outer  regions  of  the  incoming  layer  (including  the 
freestream)  curved  toward  the  surface  as  the  low-momentum  fluid  near  the  wall  is  removed  in 


Page  8 


the  spanwise  direction  by  the  main  vortical  flow.  Not  surprisingly,  the  maximum  skin-friction 
and  heat- transfer  rates  occur  near  the  jet  impingement  location. 

These  observations  are  based  almost  exclusively  on  relatively  low  Mach  number  flows  (typically 
less  than  3).  Data  at  higher  Mach  numbers,  relevant  to  hypersonic  flight,  are  virtually  non¬ 
existent,  except  for  the  recent  experimental  data  of  Bookey  et  al.  (2005)  reporting  surface  oil 
flow  visualizations  and  filtered  Rayleigh  scattering  images  for  a  Mach  8  and  Ree=3500  STBLI 
generated  by  an  8°  sharp  fin.  DNS  of  sharp  fin  interactions,  even  at  these  low  Reynolds 
numbers,  are  not  feasible,  and  robust  LES  are  required.  Currently,  LES  of  this  canonical  flow  do 
not  exist. 

1.2.  Shock  unsteadiness 

One  of  the  key  features  of  shock  wave  and  turbulent  boundary  layer  interactions  is  the  unsteady 
motion  of  the  shock.  The  shock  motion  has  a  frequency  much  lower  than  the  characteristic 
frequency  of  the  incoming  boundary  layer.  The  time  scale  of  the  low-frequency  motion  is 
0(1O5/Uoo  -  100  5/Uoo)  as  reported  in  various  experiments  such  as  Dolling  &  Or  (1985),  Selig 
(1988),  or  Dussauge,  Dupont  &  Debieve  (2006).  In  contrast,  the  characteristic  time  scale  of  the 
incoming  boundary  layer  is  0(5/Uoo). 

The  shock  unsteadiness  has  been  primarily  studied  for  two-dimensional  interactions,  where  the 
shock  translates  in  the  streamwise  direction  with  translation  magnitude  of  0(5)  and  with  smaller 
wrinkling  motion  superimposed  (Ganapathisubramani,  Clemens  &  Dolling,  2007;  Wu  &  Martin, 
2007;  Wu,  Lempert  &  Miles,  2000).  Figure  13  illustrates  such  motion.  The  spanwise  wrinkling 
is  caused  by  the  upstream  boundary  layer  structures  convecting  through  the  shock  (Erengil  & 
Dolling,  1991).  Presently,  there  are  two  schools  of  thought  that  try  to  explain  the  cause  for  the 
translation  motion,  namely  being  given  by:  (a)  the  upstream  boundary  layer,  more  recently 
'superstructures’  (Andreopoulos  &  Muck,  1987;  Erengil  &  Dolling,  1991;  Beresh,  Clemens,  & 
Dolling,  2002;  Ganapathisubramani,  Clemens  &  Dolling,  2007)  and  (b)  the  downstream 
separated  flow  (Thomas,  Putnam  &  Chu,  1994;  Dussauge,  Dupont  &  Debieve,  2006).  Recent 
analysis  on  validated  DNS  data  of  a  compression  comer  employ  correlations  of  shock  motion 
with  the  upstream  and  downstream  flow  to  find  that  the  shock  motion  is  mainly  correlated  with 
downstream  flow  dynamics  (Wu  &  Martin,  2008).  Similar  trends  are  observed  in  a  reflected 
shock  interaction  using  analyses  on  DNS  data  (Martin,  Priebe  &  Wu,  2008).  Dussauge  et  al. 
(2006)  find  that  using  StL  =  f  UU M  ,  where  L  is  the  streamwise  length  of  the  mean  separation 
bubble,  experimental  data  (covering  a  wide  range  of  Mach  numbers  and  Reynolds  numbers  and 
various  configurations)  can  be  grouped  between  StL=  0.02  and  0.05.  This  range  is  consistent 
with  the  Strouhal  number  found  in  DNS  data  (Wu  &  Martin,  2008). 

No  detailed  studies  exist  regarding  heat  transfer  or  real  gas  effects  on  shock  unsteadiness.  In 
addition,  there  exist  no  complete  and  validated  LES  data  to  study  the  long-term  dynamics  of 
shock  unsteadiness.  Recently,  Edwards,  Choi  &  Boles  (2008)  presents  a  study  of  shock  motion 
using  a  hybrid  LES/RANS  simulation  of  a  Mach  5  compression  comer  interaction.  However,  a 
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robust  and  general  methodology  for  shock-capturing  in  LES  of  STBLI  has  not  been  proven  to 
date. 

1.3.  Large  eddy  simulations 

In  the  large-eddy  simulation  (LES)  technique,  the  contribution  of  the  large,  energy-carrying 
structures  to  momentum  and  energy  transfer  is  computed  exactly,  and  the  effect  of  the  smallest 
scales  of  turbulence  is  modeled.  While  a  substantial  amount  of  research  has  been  accomplished 
for  LES  of  incompressible  flows,  applications  to  compressible  flows  have  been  significantly 
fewer.  This  is  in  part  due  to  the  limited  amount  of  accurate  and  detailed  data  at  relevant 
conditions  that  have  been  available  in  the  past  to  test  and  develop  LES  approaches  and  models. 

The  traditional  LES  approach  is  based  on  filtering  the  governing  equations  and  modeling  the 
corresponding  subgrid-scale  (SGS)  terms  (Rogallo,  1984;  Lesieur,  1996;  Piomelli,  1999; 
Meneveau  &  Katz,  2000).  The  resulting  SGS  models  are  based  on  scaling  arguments  leading  to 
eddy  viscosity  models,  (Smagorinsky,  1963;  Moin  et  al.,  1991)  physically-based  assumptions 
such  as  scale  similarity  (Bardina  et  al.,  1980),  and  practical  add-hoc  approximations  leading  to 
the  monotonically  integrated  LES  approach  (MILES)  (Boris,  1992;  Gamier  et  al,  1999;  Urbin, 
2001;  Fureby,  2002).  In  an  alternate  approach,  the  unfiltered  velocity  field  is  approximated  and 
used  to  compute  the  SGS  terms  in  their  exact  form.  Such  an  approach  is  used  in  the  estimation 
model,  where  an  estimate  of  the  unfiltered  velocity  is  obtained  by  generating  subgrid-scales  two 
times  smaller  than  the  grid  scale  through  the  nonlinear  interactions  among  the  resolved  scales 
(Domaradzki,  1997;  Domaradzki,  1999;  Dubois,  2002).  Also,  in  the  approximate  deconvolution 
model  (ADM),  a  mathematical  approximation  of  the  unfiltered  solution  is  constructed  and  used 
to  calculate  the  nonlinear  terms  in  the  filtered  governing  equations  (Stolz,  1999;  Stolz,  2001a; 
Stolz,  2001b;  von  Kaenel,  2004). 

Rizetta  and  Visbal  (2001)  perform  a  LES  of  flow  over  a  compression  comer  using  a  dynamic 
Smagorinsky  model  and  compared  their  results  to  the  DNS  of  Adams  (Adams,  2000)  and 
experimental  data  at  higher  Reynolds  numbers.  The  Reynolds  number  discrepancy  made  it 
difficult  to  draw  any  firm  conclusions  on  the  accuracy  of  the  LES.  Gamier,  Sagaut  &  Deville 
(2002)  performed  a  LES  of  a  reflected  shock  case.  Satisfactory  comparisons  in  the  mean  and 
fluctuating  streamwise  velocity  with  experiments  under  the  same  conditions  were  reported. 
However,  the  integration  time  of  the  LES  was  not  long  enough  to  study  the  low  frequency  shock 
motion.  Loginov  et  al.  (2006)  used  LES  with  ADM  to  simulate  a  compression  comer  flow 
configuration.  Comparison  with  experimental  data  showed  good  agreement  in  the  wall-pressure 
and  skin  friction  distribution.  The  comparison  of  the  velocity  profiles  in  the  streamwise 
direction  also  showed  promising  results.  However,  the  robustness  and  validity  of  the  ADM 
model  applied  to  other  interaction  configurations  and  flow  conditions  is  still  to  be  verified,  given 
the  presence  of  tunable  parameters  to  achieve  the  stability  and  accuracy  of  the  solution. 

In  the  past,  we  did  not  have  access  to  detailed  experimental  or  validated  DNS  solutions.  With 
the  refashioned  shock-capturing  methods  that  we  have  developed  in  our  group  (see  section  3), 
such  limitation  does  no  longer  exist.  We  have  access  to  true  DNS  data  of  highly  compressible 
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turbulence  interacting  with  strong  shock  waves,  and  using  these  datasets  we  can  go  back  and 
assess  the  robustness  of  numerical  methodologies  and  turbulence  models  for  LES,  making  LES 
an  indisputably  reliable  tool  for  the  prediction  STBLI. 


2.  Results  from  prior  AFOSR  (#FA9550-06-l-0323)  sponsored  research 

As  part  of  a  previous  AFOSR  grant  (#FA9550-06- 1-0323),  we  performed  a  joined  numerical  and 
experimental  effort  with  the  goal  of  validating  DNS  and  LES  for  STBLI.  In  this  regard,  detailed 
experimental  data  for  the  three  configurations  described  above  were  obtained  at  flow  conditions 
that  were  accessible  to  DNS  and  LES.  The  experiments  and  simulations  were  conducted, 
respectively,  in  the  Gas  Dynamics  Laboratory  and  the  CRoCCo  Laboratory  in  Princeton 
University.  For  the  duration  of  this  grant,  we  developed  shock-capturing  weighted-essentially 
non-oscillatory  numerical  methods  for  shock  and  turbulence  interactions,  resulting  in  robust 
shock-capturing  with  significantly  improved  grid  convergence  properties.  Currently,  shock- 
confining  filters  are  being  developed  to  enable  robust  large-eddy  simulations  of  these  flows.  The 
direct  numerical  simulation  data  match  the  experimental  data  within  the  experimental 
uncertainty.  In  turn,  the  validated  DNS  data  has  been  used  to  develop  a  greater  understanding  of 
the  upstream  and  downstream  flow  influence  on  STBLI  unsteadiness.  Below  is  a  summary  of 
results  from  previous  (past  three  years  and  a  half)  AFOSR  sponsored  research,  regarding 
numerical  methods  for  DNS,  DNS  data  quality,  validation  and  analyses,  and  large-eddy 
simulations. 


2. 1 .  Numerical  methods  for  direct  numerical  simulations 

Prior  to  this  work,  DNS  were  well-developed  for  incompressible  flows  but  not  for  highly- 
compressible  turbulent  flows.  The  simulation  of  turbulence  interacting  with  shockwaves 
requires  high  bandwidth  resolving  efficiency  and  robust  shock  capturing,  and  great  attention 
must  be  given  to  assess  the  numerical  dissipation  for  accurate  solutions.  Weighted  essentially 
non-oscillatory  (WENO)  schemes  (Jiang  &  Shu,  1996)  provide  robust  shock  capturing  for 
unsteady  flows.  With  these  schemes  the  numerical  flux  is  computed  as  a  weighted  sum  of 
candidate  stencils.  The  WENO  technique  is  best  introduced  in  the  context  of  the  one¬ 
dimensional  advection  equation,  where  the  semi-discrete  form  is 


du 

dt 


where  °  indicates  a  discrete  quantity  and  A  is  the  mesh  spacing.  The  evaluation  of  the  flux  is 
given  as  a  weighted  sum  of  candidate  stencils  by 
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where  q'k  are  the  candidate  flux  approximations  on  the  candidate  stencils  and  wh  are  the  weights 

assigned  to  each  candidate.  Figure  14  shows  the  candidate  stencils  for  an  r= 3  symmetric  WENO 
approximation.  In  perfectly  smooth  regions,  the  WENO  methodology  results  in  an  optimal 
stencil  that  can  be  optimized  for  maximum  order  of  accuracy  (Jiang  &  Shu,  1996)  or  optimal 
bandwidth  (Weirs  &  Candler,  1997).  If  a  shock  is  contained  within  a  stencil,  that  stencil  receives 
a  nearly  zero  weight,  thereby  driving  the  final  approximation  away  from  the  optimal  stencil 
while  avoiding  interpolation  across  shock  waves.  Martin  et  al.  (2006)  show  that  linearly 
bandwidth  optimized  WENO  methods  (Weirs  &  Candler,  1997),  hereby  WENOw,  can  be  used  to 
obtain  accurate  turbulence  results  for  certain  canonical  configurations.  For  isotropic  turbulence 
at  incompressible  flow  conditions,  Fig.  15a  shows  that  WENOw  schemes  can  yield  grid 
converged  results  when  using  very  fine  meshes  1283  grid  points,  relative  to  643  grid  points  for 
non-shock-capturing  schemes.  A  measure  of  numerical  dissipation  for  turbulent  boundary  layers 
is  the  skin  friction  coefficient.  For  ffeestream  Mach  numbers  in  the  range  of  3  to  8,  Fig.  15b 
shows  that  accurate  results  (within  the  uncertainty  of  the  van  Driest  II  prediction)  can  be 
obtained  with  WENOw  schemes.  Other  more  stringent  problems  could  not  be  accurately 
simulated  using  simple  linearly  bandwidth  optimized  WENO  methods.  This  is  the  case  of 
shockwave  and  turbulent  boundary  layer  calculations  or  LES  for  which  increasing  the  grid  size 
was  not  an  option.  Figures  16a  and  16b  show  the  inaccuracy  of  WENOw  for  a  STBLI  in 
comparison  with  experiments  and  a  LES  of  a  Mach  4  turbulent  boundary  layer  in  comparison  to 
DNS,  respectively.  The  inadequacy  of  simply  linearly  optimized  WENO  methods  prompted  a 
series  of  numerical  studies  to  address  the  overly  dissipative  nature  of  the  existing  WENO 
schemes. 

The  numerical  dissipation  in  WENO  comes  from  two  sources:  (1)  the  linear  dissipation  due  to 
the  theoretical  bandwidth  properties  of  the  optimal  stencil;  and  (2)  the  non-linear  dissipation  due 
to  the  non-theoretical  bandwidth  properties  of  adapted  stencils  when  deviation  from  the  optimal 
stencil  is  necessary.  The  linear  source  of  numerical  dissipation  was  addressed  by  Weirs  & 
Candler  (1997),  resulting  in  the  WENOw  schemes.  Under  AFOSR  funding  we  worked  on  robust 
solutions  to  address  the  nonlinear  numerical  dissipation  and  developed  the  linear  and  non- 
linearly  optimized  WENO  methods  (Taylor,  Wu  &  Martin,  2007;  Wu  &  Martin,  2007;  Taylor  & 
Martin,  2008).  Figure  17  plots  DNS  results  on  isotropic  turbulence  at  turbulent  Mach  number  of 
1.5  showing  that  grid  convergence  is  achieved  at  643  grid  points  using  the  new  linearly  and  non- 
linearly  optimized  WENO,  relative  to  1283  grid  points  required  using  simply  linearly  optimized 
WENO  methods.  Figure  18  shows  that  the  new  WENO  methods  result  in  accurate  DNS  data 
(that  is  within  the  experimental  uncertainty)  for  shock  wave  and  turbulent  boundary  layer 
interactions  on  affordable  grids.  We  can  obtain  converged  mean  and  fluctuating  quantities  in 
nine  days  for  the  compression  comer  shown  in  Fig.  18  using  the  new  methods,  relative  to  about  a 
year  if  were  using  WENOw.  Moreover,  the  linearly  and  nonlinearly  optimized  WENO  methods 
are  general  and  no  modifications  are  necessary  to  accurately  resolve  turbulence  or  perform  shock 
capturing.  Using  these  new  methods,  we  have  been  able  to  perform  DNS  of  STBLI,  such  as 
those  described  below,  and  compressible  isotropic  turbulence  interacting  with  imposed  shock 
waves  such  as  those  in  Fig.  2  (Taylor,  Grube  &  Martin,  2008). 
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2.2.  Validated  direct  numerical  simulation  data  of  STLBI 

Two  flow  configurations  were  considered,  namely  a  compression  comer  interaction  with 
compression  angle  of  24°  (Wu  &  Martin,  2007;  Martin  &  Wu,  2007;  Wu  &  Martin,  2008)  and  a 
reflected  shock  interaction  using  a  free  stream  flow  deflection  of  12°  (Martin,  Priebe  &  Wu, 
2008;  Priebe,  Wu  &  Martin,  2009).  For  both  cases,  the  incoming  boundary  layer  was  at  Mach 
2.9  and  Reynolds  number  based  on  momentum  thickness  of  2390,  matching  the  experimental 
data  of  Bookey  et  al.  (2005)  and  Ringuette  et  al.  (2007). 

For  the  compression  comer,  the  mean  wall-pressure  distribution  for  the  DNS  and  experimental 
data  (Bookey  et  al.,  2005)  for  the  same  configuration  and  flow  conditions  is  given  in  Fig.  19a, 
showing  good  agreement.  The  error  bars  show  an  estimated  experimental  error  of  5%.  The 
comer  is  located  at  x=0.  Figure  19b  plots  the  magnitude  of  wall-pressure  fluctuations  from  the 
DNS  data  and  experiments  (Ringuette  et  al.,  2007).  There  is  good  agreement  between  the  DNS 
and  experimental  data,  except  that  the  DNS  gives  slightly  higher  magnitude.  This  is  because  the 
synthetically  generated  turbulence  structures  in  the  initial  DNS  condition  produce  slightly  higher 
levels  of  uncorrelated  pressure  fluctuations,  or  noise,  in  the  incoming  boundary  layer.  Thus,  the 
fluctuating  wall  pressure  in  the  DNS  is  the  sum  of  the  actual  value,  p'w,  and  that  due  to 
uncorrelated  noise,  p'n  ,  and  ((pf+p^J)  ^(p'2w)  +  (p'2n) ,  since  ( 2p'wp'n)cm  be  neglected. 
An  estimate  of  the  noise  level  can  be  obtained  using  the  free  stream  value,  p\,  upstream  and 
downstream  of  the  shock  interaction  region.  The  mean  squared  of  the  pressure  fluctuations  is 
about  0.04%  and  0.16%  upstream  and  downstream  of  the  shock,  respectively.  Taking  the  square 
root  of  these  values  gives  an  amplification  factor  of  2,  and  an  rms  noise  value  of  2%  and  4% 
upstream  and  downstream  of  the  interaction,  respectively.  These  estimates  give  good 
approximations  of  the  differences  between  the  DNS  and  experimental  data  shown  in  Fig.  19b. 
The  histograms  of  the  wall-pressure  signals  for  the  DNS  and  experimental  data  at  matching 
conditions  are  shown  in  Fig.  20.  The  DNS  data  is  low-pass  filtered  at  50  kHz  to  match  the 
resolution  of  the  experiment  for  comparison.  Figure  21  plots  the  pre-multiplied  energy  spectral 
density  for  the  wall  pressure  given  by  the  same  DNS  and  experiments  in  three  streamwise 
locations:  in  the  undisturbed  boundary  layer,  at  the  mean  separation  and  at  the  first  peak  in  the 
magnitude  of  wall-pressure  fluctuation.  U„  Id  is  95  kHz  for  the  DNS  and  90  kHz  for  the 
experiments.  The  agreement  among  the  simulation  and  experimental  data  is  good,  with  the 
magnitudes  in  the  DNS  data  being  slightly  higher.  For  streamwise  locations  within  the 
separation  region  the  data  of  both  studies  show  low-frequency  peaks  at  similar  locations  0.6- 1.2 
kHz  for  DNS  and  0.6-0. 8  kHz  for  the  experiment.  These  low  peaks  correspond  to  the 
characteristic  low  frequency  of  the  shock  motion.  Both  numerical  and  experimental  spectra 
exhibit  peaks  at  high  frequencies  (of  order  105  kHz),  with  disagreement  between  the  peak 
locations  of  the  DNS  and  the  experiments.  This  is  due  to  a  combination  of  effects,  namely  the 
low-pass  filtering  of  the  experimental  signal,  which  determines  the  maximum  frequency 
resolution  at  about  17  kHz  (Ringuette  et  al.,  2007),  and  the  effect  of  the  characteristic  forcing 
frequency  imposed  by  the  rescaling  method  at  about  21  kHz  (Martin  &  Wu,  2007). 


Bookey  et  al.  (2005)  gathered  experimental  data  for  the  reflected  shock  case  for  the  same 
incoming  boundary  layer  flow  conditions  and  free  stream  deflection  angle  as  those  used  in  the 


Page  13 


DNS.  They  found  significant  three-dimensional  effects  imposed  by  the  side  walls  in  the 
experiment.  Figure  11  shows  surface  oil  visualizations  from  the  experiment  and  a  schematic 
drawing  of  the  near-wall  flow  pattern.  The  three-dimensionality  imposed  by  the  experimental 
side  walls  affects  the  flow  downstream  of  the  separation  point  and  the  wall-pressure.  The 
computer  power  required  to  simulate  the  entire  experimental  span  and  the  side  walls  renders  such 
calculations  impossible  today.  Thus,  comparing  the  DNS  and  experiment  data  is  not  sensible 
since  the  configurations  are  different.  The  numerical  method  is  general  and  robust  and  has  been 
applied  to  a  variety  of  shock  interaction  problems  over  a  range  of  conditions  without 
modification  (Wu  &  Martin,  2007;  Taylor,  Wu  &  Martin,  2007;  Taylor,  Grube  &  Martin,  2007), 
allowing  us  to  proceed  with  confidence  in  the  calculation  and  analysis  of  the  DNS  data  for  the 
reflected  shock  interaction.  Figure  12a  plots  the  wall-pressure  signals  for  the  DNS  data  at 
different  streamwise  locations,  incoming  boundary  layer,  mean  separation  point,  and  inside  the 
separated  region.  The  signals  resemble  those  for  the  compression  comer  case.  Figure  12b  plots 
the  pre-multiplied  energy  spectral  density  for  the  wall  pressure  signals.  The  peaks  associated 
with  the  characteristic  low  and  high  frequencies  are  about  0.15-0.5  kHz  and  17-40  kHz, 
respectively.  The  scaling  proposed  by  Dussauge  et  al,  StL  =  f  LI  U M ,  together  with  the  DNS 
data  for  the  compression  comer  can  be  used  to  obtain  a  theoretical  estimate  for  the  characteristic 
low-frequency  of  shock  motion  in  the  reflected  shock  case.  Figure  23  plots  the  skin  friction 
coefficient  for  the  compression  comer  and  reflected  shock  case  DNS  data,  with  x=0  at  the 
separation  point.  The  size  of  the  separation  region  for  the  reflected  shock  case  is  about  1.82 
times  that  of  the  compression  comer  case.  Using  the  scaling,  the  low  frequency  for  the  reflected 
shock  case  is  1/1.82  that  for  the  compression  comer,  or  about  0.3-0. 7  kHz,  which  is  close  to  the 
values  given  by  the  DNS. 

The  DNS  data  analysis  of  the  upstream  and  downstream  influence  on  the  shock  unsteadiness  for 
the  compression  comer  and  reflected  shock  cases  above  indicates  that  the  shock  unsteadiness  is 
driven  by  the  downstream  flow  (Martin  &  Wu,  2007;  Wu  &  Martin,  2008;  Martin,  Priebe  &  Wu, 
2008). 

2.3.  Large-eddy  simulation 

In  the  past,  we  implemented  and  validated  SGS  mixed  models  for  compressible  flows  in 
isotropic  turbulence  and  turbulent  boundary  layers  (Martin,  Piomelli  &  Candler,  2000;  Martin, 
2000)  and  a  LES  capability  in  generalized  curvilinear  coordinates  was  validated  (Martin  et  al. 
2000).  More  recently,  we  added  the  ADM  approach  to  our  LES  code  (Gmbe,  Taylor  &  Martin, 
2007).  Despite  our  LES  capability,  it  was  our  choice  not  to  perform  LES  of  STBLI  until  we  had 
validated  our  ability  to  perform  robust  DNS  of  STBLI  and  satisfactory  validation  of  the  data 
against  experiments  had  been  achieved.  Recently,  we  became  confident  that  we  had  satisfied 
this  constraint  and  we  moved  forward  in  the  development  and  assessment  of  robust  LES 
methodologies  and  models.  In  particular,  we  focused  our  efforts  in  assessing  the  need  and 
developing  shock-confining  filters  for  LES  of  turbulence  interacting  with  shock  waves  (Taylor, 
Gmbe  &  Martin,  2007;  Gmbe,  Taylor  &  Martin,  2007),  as  well  as  details  regarding  commutative 
errors  (Gmbe  &  Martin,  2009). 
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With  the  exception  of  static  eddy  viscosity  models,  most  turbulence  models  employed  in  LES 
require  explicit  application  of  filtering  operations,  sometimes  more  than  once  during  a  given  time 
advancement  stage.  Dynamic  (Germano  et  al.,  1991;  Moin  et  al.,  1991),  scale-similarity 
(Bardina  et  al.,  1980),  and  mixed  (Speziale  et  al.,  1988;  Vreman,  Geurts  &  Kuerten,  1994) 
models  rely  on  filtering  to  identify  the  smallest  resolved  length  scales;  and  ADM  relies  on  the 
iterative  application  of  filtering  to  approximately  de-filter  the  flow  solution.  Therefore  the 
calculations  of  the  unclosed  terms  of  the  filtered  Navier-Stokes  equations,  and  in  turn  the  global 
dynamics  of  the  simulated  fluid  flow,  are  directly  affected  by  the  choice  of  filtering  technique. 
In  the  vicinity  of  shocks  (or  other  discontinuities),  part  of  the  filter  stencil  may  lie  across  a  shock, 
and  in  general  this  filtering  will  cause  smearing  of  the  shock  and/or  the  creation  of  spurious 
oscillations  on  either  side  of  the  shock.  In  order  to  avoid  filtering  across  shocks,  we  have 
developed  a  shock-confining  filtering  (SCF)  technique  which  adapts  its  coefficients  in  response 
to  the  local  smoothness  of  the  flow  solution  (Grube,  Taylor  &  Martin,  2007).  Figure  24 
illustrates  schematically  the  desired  behavior  of  an  SCF,  namely  discontinuities  are  unaffected, 
but  smooth  regions  are  filtered.  In  Taylor,  Grube  &  Martin  (2007),  we  use  DNS  data  of 
shock/isotropic-turbulence  interaction  to  perform  preliminary  testing  of  the  shock-confining 
filters  against  linear  filters.  We  find  that  linear  filtering  consistently  causes  streamwise  profile 
(across  the  shock)  data  to  exhibit  anomalies  immediately  downstream  of  the  main  shock,  which 
is  indicative  of  non-negligibly  altered  global  flow  dynamics.  We  also  find  that  these  can  be 
consistently  avoided  by  the  application  of  shock-confining  filters. 

In  LES  models  using  ADM,  a  relaxation  parameter  is  introduced  in  an  attempt  to  filter  the 
spurious  oscillations  that  are  generated  during  the  linear  filtering  operation  across  a  shock  wave. 
Figure  25  shows  the  performance  of  ADM  on  laminar  shock  tube  flow  using  central  schemes  as 
in  Stolz  &  Adams  (2001).  The  generation  of  spurious  oscillations  without  the  relaxation 
parameter,  and  the  unsuccessful  removal  of  oscillations  by  the  relaxation  parameter  are  apparent. 
Figure  26a  shows  the  behavior  of  ADM  with  WENOw  not  employing  SCF,  giving  slightly  better 
results  than  ADM  with  central  schemes.  Finally,  Fig.  26b  shows  the  combination  of  linearly  and 
non-linearly  optimized  WENO  with  SCF  for  ADM,  successfully  removing  spurious  oscillations 
and  robustly  performing  shock  capturing. 

In  conclusion,  we  have  developed  shock  confining  filters  that  will  potentially  enable  robust 
(stable  without  tuning  parameters)  LES  of  STBLI.  We  have  tested  these  filters  a  priori  using 
DNS  data  of  shock  /isotropic-turbulence  interaction,  and  a  posteriori  on  decaying  compressible 
isotropic  turbulence  and  laminar  shock-tube  calculations.  Further  testing  a  posteriori  remains  to 
be  done  on  shock/isotropic-turbulence  interaction  and  STBLI.  In  addition,  further  assessment  of 
the  linear  and  non-linearly  optimized  WENO  methods  is  necessary  for  LES  resolution  in  STBLI 
to  render  such  LES  methods  general  and  robust. 
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3.  Results  from  this  AFOSR  (#FA9550-09-l-0464)  sponsored  research 

3. 1 .  Low-frequency  unsteadiness 

As  briefly  mentioned  in  Section  1.2,  a  key  feature  of  STBLI  is  strong  unsteadiness.  When  the 
flow  is  separated,  the  shock  oscillates  in  the  streamwise  direction  at  relatively  low  frequency. 
The  characteristic  frequency  of  the  shock  motion  will  typically  be  1  to  2  orders  of  magnitude 
lower,  i.e  0.1-0.01  Ux  Id  (Smits  &  Dussauge,  2006;  Dolling,  2001). 

The  cause  of  the  low-frequency  unsteadiness  is  still  under  debate.  It  has  been  proposed  that  the 
shock  motion  is  due  to  the  upstream  boundary  layer,  see  e.g.  Ganapathisubramani,  Clemens  and 
Dolling  (2007),  or,  alternatively,  that  it  is  due  to  the  downstream  separated  flow,  see  e.g. 
Dupont  et  al.  (2006),  Dussauge  et  al.  (2006),  and  Piponniau  et  al  (2009). 

Recent  large  eddy  simulations  (LESs)  and  direct  numerical  simulations  (DNSs)  capture  the  low- 
frequency  unsteadiness.  Touber  and  Sandham  (2009)  performed  the  LES  of  a  reflected  STBLI  at 
Mach  2.3  and  Re0  =5900,  matching  experimental  flow  conditions  (Dussauge  et  al.,  2006).  The 
low-frequency  unsteadiness  is  present  in  their  simulations  at  the  same  frequency  as  in 
experiments.  They  use  a  new  inflow  technique  to  ensure  that  no  artificial  low-frequency  forcing 
is  introduced  into  the  computation.  A  stability  analysis  is  also  performed,  and  this  shows  the 
presence  of  a  global  instability  mode  which  could  be  connected  to  the  observed  low-frequency 
unsteadiness.  Wu  and  Martin  (2007,  2008)  performed  the  DNS  of  a  24°  compression  ramp  at 
Mach  2.9  and  Re0  =  2300.  They  validated  their  results  (in  terms  of  separation  length,  mean  wall- 
pressure  distribution,  evolution  of  the  mean  flow  through  the  interaction)  against  experiments  by 
Bookey  et  al.  (2005a,  2005b)  at  matching  conditions.  In  addition,  Ringuette,  Wu  and  Martin 
(2008)  validated  the  fluctuating  wall-pressure  in  the  DNS  against  the  experiments  by  Ringuette 
and  Smits  (2007).  The  shock  motion  was  inferred  in  the  DNS  from  wall-pressure  fluctuations, 
and  mass-flux  fluctuations  in  the  ffeestream,  and  its  frequency  was  found  to  agree  with  a  scaling 
previously  proposed  (2006).  Priebe,  Wu  and  Martin  (2009)  performed  the  DNS  of  a  reflected 
shock  STBLI  at  Mach  2.9  and  Re0  =  2300  with  a  flow  deflection  through  the  incident  shock  of 
12°.  Both  the  simulations  by  Wu  and  Martin  (2007,  2008),  and  by  Priebe,  Wu  and  Martin  (2009) 
have  relatively  short  signal  lengths  and  cover  only  approximately  2-3  periods  of  the  low- 
frequency  unsteadiness  (based  on  its  dominant  frequency). 

During  the  Grant  effort,  direct  numerical  simulation  data  was  gathered  over  1000  6/U00  (350 
Lsep  IUX),  making  the  analysis  of  the  low-frequency  unsteadiness  possible. 
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3.2.  Flow  configuration  and  computational  set  up 


The  flow  configuration  considered  in  the  present  work  is  a  24°  compression  ramp.  The 
conditions  of  the  inflow  boundary  layer  are  Mach  2.9  and  Re0  =  2900. 

The  numerical  scheme  is  identical  to  that  used  and  validated  in  Wu  &  Martin  (2007)  in  the  DNS 
of  a  24  0  compression  ramp  flow  at  Mach  2.9  and  Re0  =  2300.  The  numerical  code  used  was 
validated  by  Wu  &  Martin  (2007)  (in  terms  of  separation  length,  mean  wall  pressure,  and 
velocity  profiles  through  the  interaction)  against  the  experiments  of  Bookey  et  al.  (2005)  at 
matching  flow  conditions.  In  addition,  the  fluctuating  wall  pressure  in  Wu  &  Martin  (2007)  was 
validated  by  Ringuette,  Wu  &  Martin  (2008)  against  the  experiments  of  Ringuette  &  Smits 
(2007).  Details  on  the  grid  size  and  resolution  appear  in  Priebe  &  Martin  (2010).  The  inflow 
condition  is  provided  by  an  auxiliary  DNS.  Details  regarding  the  boundary  and  inflow 
conditions  appear  in  Priebe  &  Martin  (2009)  and  Priebe  &  Martin  (2010),  respectively. 

3.3.  General  description  of  the  flow  field 

Figure  27  shows  a  typical  instantaneous  numerical  Schlieren  visualization  for  the  present  DNS. 
The  flow  is  from  left  to  right.  The  instantaneous  organization  of  the  flow  in  a  streamwise-wall 
normal  plane  is  apparent,  including  the  large-scale  bulges  in  the  inflow  boundary  layer,  the  shock 
wave,  and  the  strong  gradients  (large  numerical  Schlieren  index)  in  the  separated,  downstream 
flow.  A  series  of  'shocklets'  can  be  seen  in  the  flow  downstream  of  the  comer,  between  the 
boundary  layer  edge  and  the  shockwave.  These  shocklets  are  typically  associated  with  the 
turning  of  the  flow  near  reattachment. 

The  three-dimensional  instantaneous  organization  of  the  flow  is  apparent  from  figure  28,  which 
shows  an  isosurface  of  the  magnitude  of  density  gradient  (|Vp|  =  2.5p00/d),  colored  by  the 

streamwise  velocity  component.  The  color  map  is  chosen  such  that  the  color  'red'  corresponds  to 
the  ffeestream  velocity  and  'blue'  corresponds  to  reversed  flow.  It  is  apparent,  then,  that  the  flow 
at  the  comer  is  reversed.  The  large-scale  bulges  in  the  incoming  boundary  layer  are  also 
beautifully  brought  out  in  this  visualization.  Since  the  iso-surface  is  partially  translucent,  several 
shocklets  can  be  seen  in  the  flow  near  reattachment  (located  'behind'  the  shock  in  this 
visualization). 

Figure  29  shows  the  skin-friction  and  wall-pressure  distribution  along  the  streamwise  direction. 
The  abscissa  in  this  plot  is  x* ,  which  is  the  distance  along  the  wall,  rather  than  the  Cartesian 
coordinate,  v.  The  skin  friction  and  wall  pressure  shown  here  are  averaged  in  time  and  across 
the  spanwise  direction,  which  is  statistically  homogeneous.  Identifying  the  region  of  separated 
flow  with  the  region  of  negative  skin  friction,  the  separation  point  is  located  at  x*  =  -2.15, 

upstream  of  the  comer.  The  reattachment  point  is  located  at  x  =0.9 d.  The  total  separation 

,  L  =3.0(5 

length  is  thus  sep 
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The  mean  wall  pressure  distribution  displays  a  'plateau'  as  is  indicative  of  separated  flow. 

3.4.  Spectral  Analysis 

The  shock  motion  is  apparent  from  wall-pressure  signals.  Figure  30  plots  three  signals,  which 
were  obtained  at  different  streamwise  locations  from  the  DNS.  The  signals  have  been  spanwise 
averaged.  It  is  apparent  that  the  pressure  displays  a  markedly  different  behavior  at  the  three 
locations  shown.  In  the  upstream,  undisturbed  boundary  layer  (signal  shown  in  light  gray),  there 
are  no  low-frequency  oscillations,  whereas  near  the  mean-flow  separation  point  (signal  shown  in 
dark  gray),  low-frequency  oscillations  are  clearly  apparent  (and  seem  to  be  the  dominant 
contribution  to  the  r.m.s.  of  the  signal).  Further  downstream,  at  the  comer  (signal  shown  in 
black),  the  low-frequency  oscillations  still  seem  to  be  present  although  they  are  clearly  not 
dominant. 

The  length  of  the  time  history  shown  in  Figure  30  corresponds  to  the  complete  time  history 
available  from  the  DNS,  which  is  equal  to  over  10006 IUK.  When  investigating  the  low- 
frequency  unsteadiness,  the  appropriate  length  scale  seems  to  be  the  length  of  the  separated  flow 
region  (Dussauge  et  al,  2006),  L  ,  based  on  which  the  time  duration  of  the  DNS  is 

approximately  350 LsepIUx. 

Figure  31  shows  the  power  spectral  densities  of  different  wall-pressure  signals.  The  power 
spectra  are  pre-multiplied  and  normalized  by  the  rms  of  the  signals  (i.e.  the  area  under  each 
curve  is  equal  to  unity). 

Figure  31a  plots  the  power  spectrum  in  the  undisturbed  boundary  layer.  The  bulk  of  the  energy  is 
contained  in  a  broadband  peak  associated  with  the  turbulent  motions.  The  peak  is  centered  at  a 
frequency  of  approximately  UX18,  as  expected.  Superimposed  on  the  broadband  peak  are  narrow 
peaks  associated  with  the  rescaling  frequency  and  its  higher  harmonics  (the  second  and  third 
harmonic  are  visible).  These  peaks  have  compact  support  in  frequency,  and  they  contain  little 
energy  compared  to  the  turbulence.  Furthermore,  they  are  at  least  a  decade  away  from  the  low- 
frequency  range  associated  with  the  shock  motion  and  hence  do  not  pollute  the  range  of 
frequencies  that  are  of  interest  here.  We  note  that  there  is  no  energy  at  the  low  frequencies 
associated  with  the  shock  motion,  that  is,  around  values  of  St  =  fLsep  IU x  of  the  order  of  0.01. 

Figure  31b  plots  the  power  spectrum  near  the  mean- flow  separation  point.  The  broadband 
turbulence  peak  and  the  narrowband  rescaling  peaks  are  present,  but  more  than  50%  of  the 
energy  is  contained  at  much  lower  frequencies  associated  with  the  shock  motion.  The  range  of 
frequencies  found  here  agrees  with  that  determined  by  Dussauge  et  al  (2006).  They  perform  a 
survey  of  STBLI  experiments  and  find  that,  across  a  range  of  flow  conditions  and  configurations, 
the  frequency  of  the  shock  motion  lies  in  the  range  St  =  0.02  -  0.05 . 
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Figure  31c  shows  the  power  spectrum  near  the  mean- flow  reattachment  point.  Some  energy  is 
still  present  in  the  range  of  frequencies  associated  with  the  shock  motion,  but  the  bulk  of  the 
energy  has  shifted  back  to  the  turbulent  scales.  In  the  downstream  flow  (see  figure  3 Id),  the 
energy  at  the  low  frequencies  is  further  reduced,  and  almost  all  of  the  energy  is  contained  in  the 
broadband  peak  associated  with  the  turbulence  in  the  downstream  flow. 

A  different  way  of  presenting  the  power  spectra  is  given  in  figure  32,  where  contours  of  power 
spectral  density  are  shown  as  a  function  of  streamwise  distance  along  the  wall,  x* ,  and  of  non- 
dimensional  frequency.  The  contour  map  is  saturated  at  a  level  of  PSD=0.6  for  greater  clarity.  In 
the  undisturbed  boundary  layer,  the  bulk  of  the  energy  is  contained  in  the  turbulent  motions. 
Near  the  mean- flow  separation  point,  most  of  the  energy  is  contained  at  much  lower  frequencies 
associated  with  the  shock  motion.  In  the  downstream  flow,  the  bulk  of  the  energy  is  contained 
again  in  the  turbulent  scales. 

The  above  distribution  of  scales  is  in  agreement  with  the  experiments  of  Dupont  et  al.  (2006)  for 
a  low-Reynolds  number  reflected  shock  boundary  layer  interaction. 

So  far,  we  have  only  considered  signals  obtained  at  the  wall.  We  now  move  away  from  the  wall 
and  investigate  signals  inside  the  flow.  First,  let  us  look  at  the  shock  motion  in  the  ffeestream. 
We  determine  the  mean  shock  location  using  a  pressure-threshold  criterion.  At  fixed  height 
above  the  wall,  the  location  of  the  shock  is  identified  with  the  first  point  encountered  in  a 
streamwise  sweep  of  the  computational  domain  (starting  from  the  inlet),  for  which  the  pressure  is 
greater  than  1.27  pK.  Figure  33  plots  the  power  spectrum  of  the  pressure  obtained  at  the  mean 
shock  location  at  z  =  1  Ad .  Since  the  shock  moves  across  this  point,  the  bulk  of  the  energy  is 
contained  in  a  broadband  peak  in  the  range  of  low  frequencies  associated  with  the  shock  motion. 

Figure  34  shows  massflux  spectra  obtained  in  the  downstream  and  separated  flow  region.  Figure 
34a  shows  the  spectrum  of  streamwise  massflux,  obtained  at  a  location  inside  the  reversed  flow 
region  above  the  comer.  Interestingly,  the  spectrum  is  quite  'flat'.  It  consists  of  a  broadband 
peak  spanning  the  entire  range  of  investigated  frequencies.  This  indicates  that  the  motions  inside 
the  separation  bubble  have  a  large  range  of  scales.  We  note  that  there  is  a  significant  amount  of 
energy  at  the  low  frequencies  characteristic  of  the  shock  motion. 

Figure  34b  shows  the  spectrum  of  streamwise  massflux,  obtained  at  the  location  of  maximum 
rms  in  the  massflux  above  the  comer.  This  point  most  likely  lies  in  the  separated  shear  layer, 
which  lies  above  the  separation  bubble.  The  spectrum  displays  a  broadband  peak,  which  is 
centered  around  a  frequency  significantly  lower  than  the  characteristic  turbulence  frequency  in 
the  incoming  boundary  layer.  This  indicates  that  the  characteristic  frequency  of  the  turbulent 
motions  in  the  shear  layer  is  lower  than  that  of  the  energetic  eddies  in  the  incoming  boundary 
layer.  Furthermore,  the  spectra  carry  significant  energy  in  the  low  frequency  range  of  the  shock 
motion. 

3.5.  Statistical  link  between  downstream  flow  and  shock  motion 
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The  physical  mechanism  that  drives  the  low-frequency  shock  unsteadiness  is  still  under  debate.  It 
has  been  proposed  that  an  upstream  mechanism  is  at  work,  whereby  the  shock  motion  is  driven 
by  a  low-frequency  mechanism  in  the  upstream  boundary  layer.  Evidence  for  such  an  upstream 
mechanism  is  given  by,  e.g.,  Ganapathisubramani  et  al.  (2007).  In  contrast  to  this  upstream- 
driven  mechanism,  it  has  also  been  proposed  that  the  shock  motion  is  driven  by  the  downstream, 
separated  flow,  evidence  for  which  is  given  by,  e.g.,  Dupont  et  al.  (2006),  Dussauge  et  al.  (2006), 
and  Piponniau  et  al.  (2009). 


In  this  section,  we  investigate  the  statistical  link  between  the  shock  motion  and  the  upstream  and 
downstream  flow  from  the  DNS  dataset.  In  order  to  investigate  the  statistical  link,  we  consider 
the  magnitude-squared  coherence  function  (MSC).  The  MSC  between  two  signals  a  and  b  is 
defined  as 


coh(f)  = 


\PMj 

PaaWM) 


where  Paa  and  Pbb  are  the  power  spectral  densities  of  the  signals  a  and  b,  respectively,  and  Pab  is 
the  cross-power  spectral  density.  The  coherence  function  indicates,  as  a  function  of  frequency, 
how  well  signal  a  corresponds  to  signal  b.  The  MSC  function  takes  values  that  range  from  0  to  1, 
where  a  value  of  0  indicates  that  the  two  signals  are  unrelated,  and  a  value  of  1  indicates  that  a 
linear  relationship  between  the  two  signals  exists. 


As  an  indicator  for  the  shock  motion,  we  use  the  pressure  at  the  mean-shock  location  in  the 
freestream  (at  z  =  1.4(5).  The  spectrum  for  this  signal  was  shown  in  figure  33.  The  MSC 
between  the  'shock-indicator  signal'  and  several  wall-pressure  signals  in  the  separated  flow  is 
shown  in  figure  35.  As  expected,  all  plots  show  a  significant  coherence  at  the  rescaling 
frequency  and  at  its  second  and  third  harmonic.  This  behavior  is  spurious  and  unphysical,  but  we 
argue  that  it  is  acceptable  for  investigating  the  low-frequency  shock  unsteadiness,  since  the 
rescaling  frequency  lies  a  decade  above  the  dominant  frequencies  associated  with  the  shock 
motion  and  hence  does  not  pollute  the  low-frequency  part  of  the  spectra  or  coherence  plots  that  is 
of  interest  here. 


Figure  35a  shows  the  MSC  between  the  shock-indicator  signal  (signal  a)  and  the  wall  pressure 
near  the  mean-flow  separation  point  (signal  b).  There  is  no  significant  coherence  at  high 
frequencies  (except  for  the  rescaling  frequency  and  its  higher  harmonics,  as  discussed  above).  In 
the  low  frequency  range  of  the  shock  motion,  however,  there  exists  a  strong  coherence  between 
the  two  signals.  The  value  of  the  MSC  is  almost  1  in  this  range  indicating  a  near-linear 
relationship  between  the  two  signals.  We  thus  conclude  that  the  shock  motion  in  the  freestream 
and  the  wall-pressure  fluctuations  near  the  separation  point  are  related  almost  linearly,  which  is 
consistent  with  the  experimental  findings  of  Dupont  et  al.  (2006)  for  a  reflected  shock 
interaction. 

In  figures  35b,  35c,  and  35d,  the  location  of  signal  a  (the  shock-indicator  signal  in  the 
freestream)  is  kept  fixed,  while  the  location  of  signal  b  (the  wall-pressure  signal)  is  moved 
further  downstream  from  one  figure  to  the  next.  In  figure  35b,  signal  b  is  taken  under  the 
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separated  flow  region  at  a  point  downstream  of  the  shock  motion  region  but  still  upstream  of  the 
comer;  in  figure  35c,  signal  b  is  taken  near  the  mean-flow  reattachment  point;  and  in  figure  35d, 
it  is  taken  in  the  downstream  flow.  All  plots  show  strong  coherence  in  the  range  of  low 
frequencies  associated  with  the  shock  motion,  indicating  that  the  wall-pressure  under  the 
separation  bubble,  and  slightly  downstream  of  reattachment,  undergoes  low-frequency 
oscillations  that  are  coherent  with  the  shock  motion  in  the  freestream. 

3.6.  Statistical  link  between  upstream  flow  and  shock  motion 

Figure  36  shows  the  statistical  link  between  the  shock  motion  in  the  freestream  and  the  upstream 
boundary  layer.  Signal  a  is  taken  to  be  the  pressure  at  the  mean-shock  location  at  $z=1.4\delta$, 
i.e.  the  same  shock-indicator  signal  that  was  used  in  figure  35  above.  Signal  b  is  taken  to  be  the 
wall-pressure  under  the  undisturbed  boundary  layer  in  figure  36,  and  in  figure  36  it  is  taken  to  be 
the  massflux  in  the  middle  of  the  log  layer  (z+  =  60). 

In  order  to  be  able  to  draw  statistical  inferences  from  these  plots,  we  have  derived  a  99.9% 
confidence  limit,  which  is  shown  in  the  figures.  If  an  estimate  of  the  coherence  lies  above  this 
line,  then  the  probability  that  the  two  signals  are,  in  fact,  independent  (i.e.  that  the  true  coherence 
is  zero)  is  only  0.1%.  In  other  words,  if  an  estimate  lies  above  the  confidence  limit  shown,  then 
the  probability  that  the  two  signals  have  a  true  coherence  that  is  greater  than  zero  is  99.9%.  The 
confidence  limit  has  been  determined  from  Monte  Carlo  simulations. 

Both  plots  in  figure  36  show  a  weak,  but  statistically  significant  coherence  between  the  shock 
motion  and  the  upstream  flow  in  the  DNS,  in  the  sense  that  the  coherence  estimates  lie  above  the 
99.9%  confidence  limit,  at  least  for  some  frequencies  in  the  range  of  low  frequencies  associated 
with  the  shock  motion. 

3. 7.  Conclusions  from  AFOSR  Grant  #FA9550-09-l-0464 

The  direct  numerical  simulation  data  of  a  compression  ramp  shockwave  and  turbulent  boundary 
layer  covering  a  time  of  more  than  10005 IUX,  or  approximately  350 Lsep/Ux  was  gathered  to 

make  the  analysis  of  the  low  frequency  unsteadiness  possible,  and  thereby  determine  the 
upstream  and  downstream  flow  influence.  The  ramp  angle  is  24°  and  the  inflow  boundary  layer 
is  at  Mach  2.9  andRe0  =  2900. 

The  low-frequency  unsteadiness  in  the  DNS  is  apparent  from  wall-pressure  signals  near 
separation,  as  well  as  from  signals  at  the  mean  shock  location  in  the  freestream.  The  time  scales 
in  the  downstream  and  separated  flows  are  determined.  We  investigate  the  statistical  link 
between  the  shock  motion  (as  indicated  by  a  pressure  signal  in  the  freestream)  and  various 
signals  in  the  upstream  and  downstream  flow.  Strong  coherence,  at  the  low  frequencies  of  the 
shock  motion,  is  observed  with  the  downstream,  separated  flow.  These  findings  are  consistent 
with  the  experiments  of  Dupont  et  al.  (2006)  for  a  reflected  shock  interaction.  The  coherence 
with  the  upstream,  undisturbed  boundary  layer  is  shown  to  be  statistically  significant,  but  weak. 
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FIGURES 


Figure  1:  DNS  data  (Martin,  Priebe  &  Wu,  2008)  of  a  reflected  shock  interaction.  Incoming  boundary  layer  at 
Mach  2.9  and  Ree=2390,  shock  induced  by  a  12°  wedge  in  the  free  stream.  Magnitude  of  density  gradients  is 
shown. 


Figure  2.  Contours  of  normalized  density  gradient  |Vp|  /(p)  on  instantaneous  cross-sections  of  DNS  of  tested 

shock/isotropic-turbulence  interaction  configurations,  (a)  M=1.5,  Re^0=25,  Mto=0.2;  (b)  M=2.0,  Re^0=35,  Mto=0.7; 
(c)  M=2.0,  ReXo=35,  Mt0=1.3.  From  Taylor,  Grube  &  Martin  (2007). 


Page  29 


Wall  pressure,  KPa 


Figure  3:  Experimental  data  (adapted  from  Schiilein,  2006)  on  two  swept  shock  interactions  generated  by 
a  sharp  fin  on  a  Mach  5  and  Re/m=5-50xl06  boundary  layer. 


Figure  4:  Heat  transfer  as  predicted  by  various  RANS  calculations  on  the  blunt  cone  forebody  of  the 
HIFiRE-1  vehicle  (MacLean,  Wadhams  &  Holden,  2008). 
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Figure  5:  magnitude  of  species  fluctuations  in  a  reacting  Mach  4  boundary  layer  at  Ree=9480  (from  Martin  & 
Candler,  2001). 


shock 


shock 


Figure  6:  canonical  STBLI  configurations. 
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Figure  8:  Surface  pressure  distribution  on  various 
compression  corner  interactions  at  Mach  2.85 
(Settles  et  al,  1978) 


Figure  7:  Shadowgraph  images  of  Mach 
2.85  compression  corner  interactions  at 
various  wedge  angles  (Settles  et  al, 
1978) 


Figure  9:  Schematic  of  an  unseparated 

reflected  shock  interaction  (Delery  and  Marvin, 
1986). 


Figure  10:  Schematic  of  a  separated  reflected 
shock  interaction  (Delery  &  Marvin,  1986). 


Figure  11:  Three-dimensional  flow  pattern  observed  experimentally  on  a  reflected  shock  configuration:  (a)  surface 
flow  visualization  and  (b)  sketch  of  the  surface  streamline  pattern,  from  the  experiments  of  Bookey  et  al.  (2005). 
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Figure  12:  Structure  of  three-dimensional  shock-wave  boundary-layer  interaction  generated  by  a  sharp  fin  at  an  angle 
of  attack,  from  Knight  et  al.  (1987). 


0-2  0  2 


■(d) 

: 

\ 

^  * 

i1 1 

SKm 

, _  i 

•  SK? »  V  .  V 

-2 

0  2 

x/S  x/S 

Figure  13:  Contours  of  magnitude  of  pressure  gradient  showing  the  shock  location  for  two  flow  realizations  separated  by 
505/Uoo  at  distance  from  the  wall  of  25  (a,b)  and  0.95  (c,d).  Dark  indicates  large  gradient.  Straight  lines  indicate  the  location 
of  the  mean  shock  ( SKm )  and  the  instantaneous  spanwise  mean  shock  location  ( SKsm ).  From  Wu  &  Martin  (2008)  using 
DNS  data  of  a  compression  corner  interaction  with  incoming  boundary  layer  at  Mach  2.9  and  Ree=2390. 
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Figure  14:  Candidate  stencil  for  the  computation  of  the  flux  for  an  r=  3 
WENOw  scheme. 


Figure  15:  (a)  Comparison  of  DNS  data  using  WENOw  and  linear  schemes  for  incompressible  isotropic  turbulence  decay 
at  turbulent  Mach  number  of  0.1  and  Reynolds  number  based  on  Taylor  microscale  of  35;  (b)  DNS  data  using  WENOw 

schemes  for  turbulent  boundary  layers  in  the  range  of  Mach  3  to  8  and  Rer  =  pwUr8  /  =  400  in  comparison  to  van 

Driest  II  predictions  with  error  bars  at  7%. 
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Figure  16:  (a)  Comparison  between  DNS  data  using  a  WENOw  scheme  against  experiments  at  the  same  conditions  for 
a  STBLI  in  a  compression  corner  configuration.  Incoming  boundary  layer  at  Mach  3  and  Reynolds  number  based  on 
momentum  thickness  of  2400.  (b)  Comparison  between  DNS  and  LES  data  of  a  turbulent  boundary  layer  at  Mach  4 
and  Reynolds  number  based  on  momentum  thickness  of  9800  using  WENOw  schemes. 


Figure  17:  DNS  data  of  decaying  isotropic  turbulence  at  turbulent  Mach  number  of  1.5  and  Reynolds  number  based 
on  Taylor  microscale  of  50.  (a)  Density  gradients  showing  abundant  number  of  shocklets  in  the  highly  compressible 
turbulence  field;  (b)  Decay  of  turbulent  kinetic  energy  for  WENOw  and  linearly/nonlinearly  optimized  WENO 
methods  at  different  grid  resolution  showing  grid  convergence  properties. 


Figure  18:  Wall-pressure  distribution  comparing  DNS  and  experimental  data  at  the  same  conditions  and  flow 
configuration  as  in  Figure  16.  DNS  solutions  for  WENOw  and  linearly/nonlinearly  optimized  WENO 
methods.  Experimental  uncertainty  of  5%  shown.  From  Wu  &  Martin  (2007). 
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Figure  19:  Wall-pressure  distribution  for  the  compression  corner  case  with  incoming  boundary  layer  at  Mach  2.9 
and  Re0=239O  and  wall  deflection  of  24°.  (a)  Mean  from  DNS  and  experiment  (adapted  from  Wu  &  Martin,  2007) 
and  (b)  rms  from  DNS  and  experiments  (adapted  from  Ringuette,  Wu  &  Martin,  2008). 
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Figure  20:  Wall-pressure  signal  from  DNS  data  at  three  streamwise  locations  compared  to  the  experiments  at 
the  same  flow  conditions  as  Fig  19.  (a)  DNS  signal.  The  DNS  data  are  low-pass  filtered  at  the  same  cut-off 
frequency  as  the  experiments,  (b)  Experimental  data  at  matching  conditions.  Dashed  lines  indicate  mean  values. 
Adapted  from  Ringuette,  Wu  &  Martin  (2008). 


Page  37 


Figure  21:  Pre-multiplied  energy  spectral  density  of  the  wall  signal  at  three  different  streamwise 
locations  for  the  DNS  (lines)  and  experiments  (symbols)  at  the  same  flow  conditions  as  in  Fig  19. 
The  streamwise  locations  correspond  to  the  incoming  boundary  layer,  the  mean  separation  point, 

and  the  peak  in  the  p'wRMS  curve,  respectively.  Adapted  from  Ringuette,  Wu  &  Martin  (2008). 
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Figure  22:  DNS  data  of  a  reflected  shock  configuration  with  incoming  boundary  layer  at  Mach  2.9  and 
Re0=239O  generated  by  a  12°  wedge  in  the  freestream.  (a)  Wall-presure  signal  and  (b)  pre-multiplied 
energy  spectral  density  for  the  wall-pressure  at  different  streamwise  locations.  Adapted  from  Martin,  Priebe 
&  Wu  (2008). 
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Figure  23:  Skin  friction  coefficient  from  DNS.  Adapted  from  Martin,  Priebe  &  Wu  (2008). 
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Figure  24:  Ideal  behavior  of  a  shock-confming  filter. 
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Figure  25:  Density  profiles  in  shocktube  problem  from  WENO-based  DNS,  ADM  with  central  differencing,  and 
ADM  with  central  differencing  and  regularization.  From  Grube,  Taylor  &  Martin  (2007). 
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Figure  26a:  Density  profiles  in  shocktube  problem  from  WENO-based  DNS,  and  WENO-based  ADM  using  a 
linear  tophat  filter.  From  Grube,  Taylor  &  Martin  (2007). 
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Figure  26b:  Density  profiles  in  shocktube  problem  from  WENO-based  DNS,  and  linearly  and  non-linearly 
optimized  WENO  ADM  using  shock-confining  tophat  filter.  From  Grube,  Taylor  &  Martin  (2007). 
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Figure  28:  Typical  instantaneous  numerical  Schlieren 
visualization  of  a  streamwise-wall  normal  plane 


Figure  29:  Instantaneous  flow  visualization. 
Isosurface  of  density  gradient 

(|v p\  =  2.5 / 8^ ,  colored  by  the  streamwise 

velocity  component ,  where  the  color  red 
corresponds  to  the  freestream  velocity,  and  blue 
corresponds  to  reversed  flow. 
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Figure  29:  Mean  distribution  of  skin  friction  (top)  and  wall  pressure  (bottom) 
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Figure  32:  Contour  map  of  premultiplied 

power  spectrum.  The  variable  X  is  the 
streamwise  distance  from  the  comer, 
measured  along  the  wall. 


st  =  fLsep/u. 


Figure  33:  Pre-multiplied  spectmm  of  pressure  at 
mean  shock  location  in  the  freestream 


(z  =  1.46). 
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Figure  34:  Pre-multiplied  spectra  of  streamwise  massflux  in  (left)  the  reversed  flow  region  above  the  corner 
and  (right)  the  shear  layer  above  the  separation  bubble. 
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Figure  35:  Statistical  link  between  the  shock  motion  and  the  separated  flow.  MSC  between  the  pressure  at  the 
mean  shock  location  in  the  freestream  ($z=1.4{\delta}$)  and  the  wall  pressure  at  different  streamwise 
locations:  (a)  near  the  mean-flow  separation  point,  (b)  in  the  separated  flow  region  upstream  of  the  comer,  (c) 
near  the  mean-flow  reattachment  point,  and  (d)  in  the  downstream  flow. 
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Figure  36:  Statistical  link  between  the  shock  motion  and  the  upstream  flow.  MSC 
between  the  pressure  at  the  mean  shock  location  in  the  freestream  (z  =  1 .4(5)  and  (top) 
the  upstream  wall  pressure,  and  (bottom)  the  massflux  in  the  middle  of  the  log  layer 
^Z+  =  60  jin  the  upstream  boundary  layer.  See  text  for  a  derivation  and  explanation  of  the 
confidence  limit  shown  in  the  figures. 
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